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Abstract 

Background: About 60% of Pheochromocytoma (PCC) and Paraganglioma (PGL) patients have either germline or somatic 
mutations in one of the 12 proposed disease causing genes; SDHA, SDHB, SDHC, SDHD, SDHAF2, VHL, EPAS1, RET, NF1, 
TMEM127, MAX and H-RAS. Selective screening for germline mutations is routinely performed in clinical management of 
these diseases. Testing for somatic alterations is not performed on a regular basis because of limitations in interpreting the 
results. 

Aim:lhe purpose of the study was to investigate genetic events and phenotype correlations in a large cohort of PCC and 
PGL tumours. 

Methods: A total of 101 tumours from 89 patients with PCC and PGL were re-sequenced for a panel of 10 disease causing 
genes using automated Sanger sequencing. Selected samples were analysed with Multiplex Ligation-dependent Probe 
Amplification and/or SNParray. 

Results: Pathogenic genetic variants were found in tumours from 33 individual patients (37%), 14 (16%) were discovered in 
constitutional DNA and 16 (18%) were confirmed as somatic. Loss of heterozygosity (LOH) was observed in 1/1 SDHB, 11/11 
VHL and 3/3 NF1 -associated tumours. In patients with somatic mutations there were no recurrences in contrast to carriers of 
germline mutations (P = 0.022). SDHx/VHUEPASl associated cases had higher norepinephrine output (P = 0.03) and lower 
epinephrine output (P<0.001) compared to RET/NF1 /H-RAS cases. 

Conclusion: Somatic mutations are frequent events in PCC and PGL tumours. Tumour genotype may be further investigated 
as prognostic factors in these diseases. Growing evidence suggest that analysis of tumour DNA could have an impact on the 
management of these patients. 
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Introduction 

Pheochromocytoma (PCC) and paraganglioma (PGL) are rare 
neural crest-derived tumours arising in the adrenal medulla (PCC) 
or autonomic ganglia (PGL). A majority of patients present with a 
focal tumour lesion and may be cured with RO resection [1-3]. 
However, even following an apparently successful surgical 
resection, there may be a risk of local or metastatic recurrence 
that motivates a long follow up period [2,4,5]. Age, familial disease 
and tumour size correlate to increased risk of malignancy and 
recurrence [2] and histologic criteria may also aid in predicting 
risk for malignant disease [6,7]. Translational studies show that 
approximately 60% of PCC and PGL cases have either germline 
or somatic mutations in one of 13 suggested disease causing loci; 



SDH subunits A, B, C and D, SDHAF2, VHL, EPAS1, RET, MF1, 
TMEM127, MAX and H-RAS [8-20]. In the clinical setting, 
genetic screening of these genes by fragment prioritization of 
germline DNA is regarded as golden standard of care, and may 
have a substantial impact on patient management [21,22]. 
Depending on the affected gene, the risk of local recurrence 
and/or metastatic disease can be estimated and guide in the 
selection of appropriate preventive measures [22]. Screening for 
tumour specific genetic events is not recommended in clinical 
practice due to a lack of genotype-phenotype correlations that 
could justify the necessary resource allocation [8,23]. However, as 
most genetic screening studies of PCC and PGL tumours have 
been performed on small cohorts, we hypothesized that a 
comprehensive genetic screening in a large clinically annotated 
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cohort could be informative. The aim of this study was, to describe 
the genetic landscape of PCC and PGL tumours and to correlate 
tumour genotype with patient characteristics and disease outcome. 

Patients and Methods 

Patients 

This is a single centre, retrospective study of 101 tumour 
samples from 89 patients with PCC and PGL treated at the 
Department of Surgery, Uppsala university hospital, Sweden. 
Selected patients were previously screened for mutations in H-RAS 
describing somatic genetic variants [n = 4) and MAX describing no 
pathogenic genetic variant [20,24]. Fourteen patients were 
clinically diagnosed with hereditary syndromes; familial paragan- 
glioma type 4 (PGL4; n = 2), Von Hippel Lindau syndrome (VHL; 
n = 4) and Multiple Endocrine Neoplasia type 2 (MEN2; n = 8) by 
certified genetic testing laboratories. Four additional patients had 
been diagnosed with Neurofibromatosis type 1 (NF1) by presence 
of clinical criteria. DNA samples from 195 healthy and unrelated 
individuals were utilized as a control to determine frequency of 
germline variants with unknown significance (VUS) in Swedish 
population. 

Ethical Statement 

Ethical approval was obtained from the regional ethics committee 
in Uppsala as well as written informed consent from the individual 
patients. All patients were above 18 years of age at the time of 
inclusion. 

Clinical Data 

Age at diagnosis was set at the time for radiological diagnosis. 
Tumour size was calculated as the mean of two diameters. Cases 
were classified as metastatic by the presence of invasion into 
nonchromaffin organs determined by histological examination or 
radiological/molecular imaging. Preoperative urinary catechol- 
amines measurements analysed in a clinical setting were included 
for evaluation. Norepinephrine assays had been performed with two 
different reference intervals, <350 nmol/24 h and <400 nmol/ 
24 h depending on utilized assays. Reference intervals for urinary 
epinephrine was <90 nmol/24 h, for plasma normetanephrine < 
0,6 nmol/L and for plasma metanephrine <0,3 nmol/L. 

DNA Extraction & Sequencing 

DNA was extracted from cryosections of tumour samples, 
peripheral blood and/or normal tissue, using DNeasy Blood & 
Tissue Kit (Qiagen, Hilden, Germany) as previously described [25] . 
Cryosections from included tumours were analysed for tissue 
morphology and selected samples were macrodissected in order to 
reduce contamination of normal cells. Using a phenotype guided 
fragment prioritization approach [1,26], exons and intron-exon 
boundaries of SDHB (NMJ303000.2), SDHC (NMJJ03001.3), 
SDHD (NM_003002.2), SDHAF2 (NMJD17841.2), VHL (NM_ 
000551.3), EPAS1 (exons 9 and 12, NM_00 1430.4), RET (exons 
10-11 and 13-16, NM_020975.4), TMEM127 (NM_0 17849.3), 
MAX (NM_002382.3) and H-RAS (exons 2 and 3, NM_1 76795.3) 
were amplified by PCR and sequenced using automated Sanger 
sequencing (Beckman Coulter Genomics, Takeley, UK). In patients 
with pathogenic germline variants, all exons and intron-exon 
boundaries of the affected gene were sequenced in order to 
investigate a potentially inactivating variant on the 2 nd allele. 
Ultimately, patients without detected pathogenic mutations in this 
study had been screened for all ten above mentioned genes. Primer 
sequences can be obtained by request. 



Mutational Analysis 

Chromatograms generated by Sanger Sequencing were reviewed 
using CLC genomics workbench 5.5 (CLC bio, Aarhus, Denmark). 
Genetic variants were annotated for overlapping information 
available in public databases; Catalogue of Somatic Mutations in 
Cancer (COSMIC) [27], the Single Nucleotide Polymorphism 
database (dbSNP), Human Genome Mutation Database (HGMD 
public) [28] and Leiden Open source Variation Database (LOVD). 
In silico analysis was performed using Sorting Intolerant From 
Tolerant (SIFT) [29] and Polymorphism Phenotyping v2 (Poly- 
phen-2) [30]. 

Multiplex Ligation-dependent Probe Amplification 

Inclusion criteria for analysis were absence of a pathogenic 
germline variant. Only DNA extracted from blood/ normal tissue 
was selected for Multiplex ligation-dependent probe amplification 
(MLPA). Included samples were analysed with SALSA MLPA 
P226 SDH and P016-B2 VHL probe mixes (MRC-Holland, 
Amsterdam, Netherlands) that have coverage of the SDHA, SDHB, 
SDHC, SDHD, SDHAF2 and VHL loci. Reactions were carried out 
as previously described [3 1] and an ABI3 1 30x1 Genetic Analyzer 
(Life Technologies, Carlsbad, CA) was used for fragment 
separation. The MLPA data were analyzed using GeneMapper 
4.0 genotyping software (Applied Biosystem) and SeqPilot version 
3.3.2 (JSI medical systems GmBH, Kippenheim, Germany). The 
experiments were carried out at a laboratory certified for clinical 
use and analysed by an experienced clinical investigator (MN). 

Single Nucleotide Polymorphism Array 

Inclusion criteria were presence of pathogenic or unknown 
variants in SDHB, SDHC, VHL or clinical criteria of NF 1 . Tumour 
DNA from the selected samples were subjected to SNParray 
analysis; using Illumina Omnil-Quad or Omni2,5-Quad chips 
(Illumina Inc, CA, USA), containing 1,140,419 and 2,379,855 
probes respectively. Hybridization and sequencing was performed 
by university core facilities, SNP&SEQ, Technology Platform in 
Uppsala, Sweden (http://molmed.medsci.uu.se/SNP+SEQ+Technology 
+Platform/). Generated data was imported into the Illumina BeadStudio 
(Illumina, CA USA) software and analysed using Nexus Copy Number 
Variation 7.0 Build 7887 (Biodiscovery Inc, CA USA) to detect copy 
number variation and allelic imbalance using default thresholds The faction 
of DNA derived from tumour cells were estimated by analysing the B-allele 
frequency in regions showing copy number alterations [32]. Samples with 
tumour cell purity <70% were analysed with adjusted settings. 

Statistics 

SPSS 19 (IBM, Armonk NY, US) was used for statistical 
calculations. Chi 2 test was performed for analysis of nominal 
variables. As age at diagnosis, tumour size and catecholamine 
output were not normally distributed, hence a non-parametric test 
(Mann- Whitney U test) was selected for analysis of these scaled 
values. Patients with variants of unknown significance (VUS) were 
classified as wild type. Multiple regression test of phenotype 
correlation to carrier status and genotypes were not possible to 
perform due to the low number of observations. P-values <0.05 
were considered as significant and <0.1 as borderline significant. 

Results 

Patient characteristics are described in table 1 . The median age 
at diagnosis was 49 years (range 15-85) and there were 35 males 
and 53 females. Eighty patients had adrenal PCC (31 left adrenal, 
37 right adrenal, 10 bilateral), eight had thoracoabdominal PGLs 
and there were one head and neck PGL. The median tumour size 
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was 55 mm (range 8-170). There were wide discrepancies in 
catecholamine output; urinary norepinephrine range 1 1 2- 
19130 nmol/24 h (ref <400/<350 nmol/24 h), urinary epineph- 
rine range 19-33322 nmol/24 h (ref <90 nmol/24 h), plasma 
normetanephrine range 1-37 nmol/L (ref <0,6 nmol/L), and 
plasma metanephrine range 0-140 nmol/L (ref <0,3 nmol/L). 
There were 13 patients with recurrent disease, eight of whom 
recurred with distant metastases and five with local recurrences. 
One additional patient had distant metastases at the time of 
diagnosis. The median follow up time was 106 months (range 0- 
714 months). 

Genetic Screening 

A total of 33 patients (37%) had a pathogenic genetic variant in 
included disease causing loci. We could confirm 14 of these 
mutations in constitutional DNA (Figure 1 a) and 1 6 were confirmed 
as somatic by absence in constitutional DNA (Figures lb and c). 
There was no constitutional DNA available for three of the patients 
with pathogenic mutations. All genetic variants previously clinically 
diagnosed by the Department of Clinical Genetics (re = 14) were 
verified in the present study. There were four additional patients 
with clinical criteria of Neurofibromatosis type 1 but with no genetic 
testing. Cases with discovered mutations and corresponding clinical 
characteristics are presented in Table 2. 

Two related patients (mother and son) presented with multiple 
abdominal PGL. Both had several local recurrences that were not 
classified as metastatic lesions. Re-sequencing revealed a patho- 
genic nonsense mutation in SDHB; c.268C>T, p.Arg90* [13] that 
was present in DNA from blood. 

There were 12 cases with pathogenic mutations in the VHL 
gene. Patient number 9 with bilateral PCC (index case) had a 
pathogenic germline missense mutation in VHL; c.482G>A 
p.Argl61Gln. A family comprising of three siblings with bilateral 
or unilateral PCC had pathogenic germline missense mutation in 
VHL; c.499C>T, p.Argl67Trp [33]. Both p.Argl61Gln and 
p.Argl67Trp had previously been reported as pathogenic [33]. 
Seven patients had somatic mutations in VHL. There were six 
unique SNVs; c.193T>G, p.Ser65Ala; c.l93T>A, p.Ser65Thr; 
c.238A>G, p.SerSOGly; c.458T>A, p.Leul53Gln; c.475A>G, 
p.Lysl59Glu and c.551T>A, p.Leul84His in one patient each. 
Patient number 4 had a 30 base pair deletion c. 163_192del, 
p.Glu55_Arg64del that was absent in DNA from peripheral blood. 
All carriers of somatic VHL mutations had unilateral PCC, 
sporadic disease presentation and there were no apparent signs or 
symptoms of VHL syndrome. 

Two patients had mutations in EPAS1 which were absent in 
DNA from their blood; one c.l586T>C, p.Leu529Pro and one 
C.15890T, p.Ala530Val. These mutations are previously 
described as pathogenic [16,34]. Both patients had sporadic 
disease presentation. Patient 1 6 had borderline polycytemia with a 
haemoglobin level of 150 g/L (reference interval 120-150 g/L). 

There were 13 patients with pathogenic mutations in RET. 
Eight had germline pathogenic mutations and clinical character- 
istics of MEN2 syndrome; c.1826G>C, p.Cys609Ser; c.1832G> 
A, p.Cys611Tyr; c.l900T>C, p.Cys634Arg; c.1900T>G, 
p.Cys634Gly; c.1901G>A, p. Cys634Tyr; c.2410G>A, p.Val804- 
Met in one patient each and c.2753T>C, p.Met918Thr in tow 
different patients. Two patients with unilateral PCC and sporadic 
disease presentation had somatic mutation in RET; c. 1900T>G, 
p.Cys634Gly and C.1900T>C, p.Cys634Arg. For three of the 
patients with SNVs in RET there were no constitutional DNA 
available; c.l891G>T, p.Asp631Tyr in one patient and 
c.2753T>C, p.Met918Thr in two patients. All three cases had 
sporadic disease presentation and there were no signs or symptoms 



suggesting MEN2 syndrome. All these RET mutations are 
described as pathogenic in the literature [8,33]. 

Four patients had previously been described with somatic H- 
RAS mutations [20]. One additional somatic mutation in H-RAS; 
c. 181 OA, p.Gln61Lys; was detected in a male patient that had 
sporadic disease presentation. 

Variants of Unknown Significance 

Four patients had germline variants of unknown significance 
(VUS). There were two VUS in SDHC; c.328C>T, ProllOSer 
[35] and c.490A>T, Metl64Leu in two different patients with 
unilateral PCC and sporadic disease presentation. Both mutations 
were available in constitutional DNA. Succinate dehydrogenase submit 
C Metl64Leu has been reported to have impact in functional 
models but was classified as benign in vivo [36]. The pathogenicity 
of ProllOSer has not been investigated in detail [35]. Succinate 
dehydrogenase submit C Codon 110 is conserved among vertebrate 
orthologs whereas codon 164 is not a conserved residue. In silico 
analysis determined the variants as benign, SIFT (0,93 and 0,96) 
and Polyphen2 (0,231 and 0,0). 

Patient 36, a 28 year old male, was referred to the local hospital 
due to hypertension, headache and multiple episodes of syncope. 
Urinary norepinephrine was elevated, 1752 nmol/24 h (<350), 
but epinephrine was within reference margins 37 nmol/24 h (< 
90). A computed tomography revealed an abdominal mass located 
in close proximity to the left renal artery and vein. The lesion was 
surgically resected and the pathology report showed a PGL, 
20 mm in size with a Ki67 index of <1%. Genetic screening of 
RET as well as Succinate Dehydrogenase subunits B and D was normal. 
The proband remains free of recurrence 2 1 months following the 
initial diagnosis. There were no apparent signs or symptoms of von 
Hippel Lindau syndrome. Sanger sequencing revealed a single 
nucleotide polymorphism in VHL; c.548C>T, p.Serl83Leu. 
Multiplex ligation-dependent probe amplification of constitutional 
DNA did not show any pathological imbalances. In tumour tissue, 
loss of heterozogosity and copy number loss was observed on the 
whole arm of chromosome 3p by Omni- 1 -quad SNP array 
(Ulumina, CA, VSA).VHL; c.548C>T, p.Serl83Leu has been 
classified as pathogenic in a functional model but the individual 
contribution of the allele to patient phenotype is not fully described 
[37]. Screening of 190 healthy individuals revealed homozygous C 
allele in all cases. Codon 183 is conserved among mammalian 
orthologs and in silico analysis determined the variant as probably 
pathogenic: SIFT (0,18) and Polyphen2 (1,0). 

Patient 37 was diagnosed with a unilateral PCC at 27 years of age 
having apparently sporadic presentation. Re-sequencing revealed a 
RET mutation c.2372A>T, p.Tyr791Phe that was found in 
constitutional DNA. The pathogenicity of RET p.Tyr791Phe is 
disputed [38,39]. 

There were no pathogenic variants discovered in SDHAF2, 
TMEM127 and MAX. 

Structural Variations 

Multiplex ligation-dependent probe amplification analysis did 
not reveal any copy number gains or losses. Analysis of SNP array 
data showed loss of heterozygosity (LOH) located to the genomic 
coordinates of the respective gene in 1 / 1 SDHB, 11/11 VHL and 
3/3 NFl-related tumours (table 3, figure 2). There were no LOH 
at coordinates corresponding to SDHC loci in tumours from 
patients with germline SDHC Prol lOSer and Metl64Leu variants. 
Patient 36 with a germline VHL />.Serl83Leu had LOH at the 
VHL loci. Analysis failed due to corrupted data in two cases, one 
SDHB and one jVFi-related tumour. 
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Figures 1. Chromatograms exported from CLC Genomics Workbench 5.5 displaying (A) Pathogenic genetic variants available in 
constitutional DNA, (B) confirmed somatic variants and (C) 30 base pair somatic deletion in VHL 

doi:1 0.1 371 /journal.pone.0086756.g001 



Statistical Correlation 

Statistical analyses of carrier status and genotype correlations to 
phenotype are presented in Tables 1 and 4. Stratified into groups 
accordingly to mutation status, germline carriers had a age at 
diagnosis that were significantly lower (median 29,5 years) 
compared to those with somatic aberrations (median 48 years, 
P — 0.025) as well as patients without known mutations (median 53 
years, P— 0.002). The frequency of mutifocal tumours were also 
different in germline carriers (53%) compared to patients with 
somatic carrier status (0%, P<0.001) as well as those without 
known mutations (2%, P<0.001). 

No cases of recurrent disease were observed in somatic carriers 
(0%), in contrast to germline carriers (28%, P = 0.022) and 
borderline significant compared to patients without discovered 
mutations (15%, P = 0.07). Preoperative levels of urine norepi- 
nephrine were lower in patients with germline carrier status 
compared to somatic carriers and those without mutation 
(P— 0.049 and P— 0.033 respectively). Gender, tumour size, 
tumour localization, metastatic disease as well as urine and plasma 
epinephrine output were not different among the three carrier 
status groups. 

Stratification according to genotype into cluster 1 ; SDHx/ VHL/ 
EPAS1 mutants and cluster 2; RET/ NF1 / H-RAS mutants, resulted 
in a difference in age at diagnosis between cluster 2 carriers 
(median 45) and patients without mutations (median 53, 
P — 0.036). A borderline significance difference was also noted 
for PCC localization, there were one left adrenal and eight right 
adrenals affected in cluster 1 compared to eight left and eight 
right-sided tumours in cluster 2 (P — 0.052). A difference in PCC 
lateralization was also noted between cluster 1 patients and those 
without discovered mutations (P— 0.028). In patients without 
mutation there were one case with multiple PGL and none 
bilateral PCC, different than in cluster 1 (4 cases total, P<0.001) 
and cluster 2 (5 cases total, P — 0.004). Urine norepinephrine levels 
(Figures 3a and b) were higher in cluster 1 patients (median 2439/ 
24 h, P=0.03) and those without mutation (median 2181,5, 
P — 0.018) compared to cluster 2 (median 862). Reversely 
Epinephrine levels were lower in cluster 1 (median 58/24 h, P< 
0.001) and in those without mutation (median 247 nmol/24 h, 
P— 0.002) compared to cluster 2 patients (median 520 nmol/ 
24 h). Age at diagnosis, gender, plasma cathecolamines as well as 
recurrent and metastatic disease were not different among the 
three groups. 



Discussion 

Genetics 

We have analysed 101 PCC and PGL tumours for SNVs in nine 
different genes, complemented by selective MLPA and SNP array 
analysis. A total of 33 patients (37%) had pathogenic variants in 
SDHB, VHL, EPAS1, RET and H-RAS. Including patients with 
clinical criteria of Neurofibromatosis type 1 and loss of heterozy- 
gosity at the NF1 locus, 41 % of the cohort could be associated with 
genetic aberrations in known genes. We did not find any 
pathogenic germline mutations that had not been previously 
discovered by clinical screening, this low frequency of germline 
mutations in apparently sporadic patients have previously been 
described in Swedish patients [40]. Considering the characteristics 
of the cohort with a predominance of benign and unilateral PCC 
having sporadic presentation, the frequencies of pathogenic 
genetic variants and patient carrier status were similar to those 
previously reported [8]. Loss of heterozygosity could be detected in 
tumour DNA from 11/12 patients having somatic or germline 
mutations in the SDHB or VHL genes. SNParray analysis of 
tumour DNA from patient 10 (VHL p.Argf61Gln) did not show 
LOH at any locus. This tumour sample will have to be carefully 
reviewed for contaminating wild type cells and re-analysed. 
Patients with clinical criteria of NF1 syndrome did not have a 
clinical diagnostic genetic test performed, probably due to the size 
of the NF1 gene and the number of possible loci. Three of these 
patients tumours were analysed with SNParray that revealed loss 
of heterozygosity at the NF1 loci in tumour DNA in all investigated 
cases. This strongly suggest that these patients do have a germline 
mutation in MF1 [15,18]. The interpretation of clinical correla- 
tions in sporadic patients presented by this study is limited by the 
absence of analysis of the MF1 gene that was recently found to be 
commonly affected in patients with sporadic PCC and PGL 
[15,18]. To analyse this extensive loci, further studies may utilize 
Next Generation Sequencing or SNP arrays [35,41]. 

Multiplex ligation-dependent probe amplification analysis of 
constitutional DNA showed no copy number gains or losses in any 
of the investigated cases. Both the MLPA laboratory workflow and 
result analysis were performed using robust workflows by 
experienced clinical investigators, ensuring high reliability of these 
results. 
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Variants of Unknown Significance 

Multiple factors contribute to determine the disease causing 
impact of a specific genetic variant; variant deleteriousness, 
probands phenotype, family history, molecular and in silico 
characterization as wells as the allele frequency in a population 
without disease [42]. For SDHC ProllOSer and Metl64Leu, 
clinical presentation and family history did not indicate familial 
paraganglioma type 3. Analysis of the biochemical phenotype 
revealed a high norepinephrine to epinephrine ratio, pointing to a 
disease causing mutation in the SDHx or VHL loci [43] . Available 
literature did not support neither classification of pathogenic nor 
benign status of SDHC ProllOSer and Metl64Leu. In silico 
analysis determined these variants to have a benign effect on 
protein structure. But as in silico methods have a high rate of false 
negative predictions [44] and the fact that the alleles are not 
reported in any normal population we classified the variants as 
VUSs. 

In patient 37 with the germline VHL p.Serl83Leu variant there 
was no family history suggesting VHL syndrome. The patient 
phenotype with a focal PGL did not strongly indicate a germline 
mutation even though up to 10% of patients may present with 
PGL [26] . The ratio of norepinephrine to epinephrine was high 
indicating a mutation in SDHx or VHL genes [43]. Screening of 
190 healthy subjects did not find this variant nor was it found by a 
database search, supporting that the variant is not a common 
polymorphism. In silico calculation determined the variant as 
probable pathogenic. This variant's contribution to patient disease 
was classified as unknown. However, several indicators point out 
this variant as potentially pathogenic and as most variants in the 
VHL gene are pathogenic, this variant should influence the clinical 
management of the proband. 

Genotype Clustering 

The molecular phenotype of EPAS1 mutated tumours has been 
disputed [16,17]. However, clinical correlations and experimental 
investigations have indicated similarities of EPAS1 lesions to SDHx 
and VHL mutated tumours. Catecholamine production detected in 
EPAS1 carriers in this study suggested cluster 1 differentiation. The 
available literature strongly suggests a common pathway for H- 
RAS with cluster 2 genes RET and NF1 [20,43]. Additionally, the 
catecholamine output observed in H-RAS mutated tumours 
resemble that of RET and MF1 [43]. We included EPAS1 in 
cluster 1 and H-RAS in cluster 2 as well as patients with clinical 
criteria of NF1. 

Genotype Phenotype Correlations 

As previously described, germline carriers were significantly 
younger at time of diagnosis, and had a higher frequency of 
multifocal disease [26]. Among carriers with somatic mutations 
there were no recurrences, nor were there any cases of metastatic 
disease. Stratified into cluster 1 and 2 genotypes, differences in 
biochemistry output was observed, in line with previous studies 
investigating germline carriers [43]. A difference in the location of 
adrenal tumours was noted with cluster 1 patients having the right 
adrenal affected in all but one case. No obvious explanation to this 
lateralisation was found in the literature and studies of larger 
cohorts is needed to confirm this observation. 

Five of 13 patients with recurrent disease were carriers of 
germline mutations in SDHB or RET, the remaining eight patients 
had sporadic disease presentation. The frequency of discovered 
variants in included loci were different in metastatic (22%) 
compared to non-metastatic cases (40%). Further genetic investi- 
gations of sporadic patients with recurrent and/or metastatic 
tumours could help to identify novel disease causing genes 
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Table 3. Copy number variation and loss of heterozygosity by SNP array. 
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Tumours with loss of heterozygosity at loci of mutated tumour suppressor. Patients with diagnostic criteria of NF1 were considered as potential carriers of a pathogenic 
variant in the NF1 gene. NA, Not Available. 
doi:1 0.1 371 /journal.pone.0086756.t003 



associated with a high risk of recurrence. Genetic prognostic 
markers could potentially be of value in a clinical context in order 
to personalize the extent of follow up [5,45]. Planned studies 
aiming at identifying factors that might predict recurrent disease 
would favourably include tumour genotype characterization. 

Tumour Genotype in Clinical Management 

There is a growing rationale for analysing somatic events in 
PCC and PGL tumours as a diagnostic test: (1) EPAS1 mutations 
may occur early in embryogenesis and these mosaic carriers are 
not found by analysis of DNA in peripheral blood [16,46]; (2) 
Translational studies have suggested using genotype as predictive 



markers for sensitivity to targeted therapy; SDHx/ VHL mutations 
might benefit from antiangiogenic treatment whereas RET/ MF1 / 
TMEM127/MAX driven tumours could benefit from inhibitions of 
kinase pathways [47,48]. The present study further suggests that 
tumour genotype might have additional prognostic implications. 
Even though methods for using formalin fixed archived tissue are 
evolving rapidly, analysis of tumour DNA is favourably performed 
using high quality fresh frozen genetic material thus discussions 
regarding the current routines, archiving only formalin fixed 
tissue, are warranted. 



3p26.3-q29 17p13.3-q25.3 

VHL NF1 

I I 

OCniMlZCIin) (TUTSI) .Allelic imbalance 

* Variant of unknown significance 



Copy number loss 
Copy number gain 
LOH 



Figure 2. Detected copy number events from SNP array as displayed by Nexus copy number 7.0. Results are separated by chromosome 
and presented for each individual tumour with patient id at the left margin. Presented data constitute cases harbouring pathogenic genetic variants 
in VHL (n = 1 1) as well as patients with clinical criteria of NF1 (n = 3). Data is also presented for patient 36 that harboured a VHL mutation of unknown 
significance. Colour annotation indicates copy number loss (red), copy number gain (blue), loss of heterozygosity (yellow) and allelic imbalance 
(magenta). Arrows indicates VHL (chromosome 3) and NF1 (chromosome 17) loci. Loss of heterozygosity at VHL locus was detected in 1 1/1 1 tumours 
with pathogenic VHL mutations and in 3/3 tumours from patients with clinical criteria of NF1. 
doi:1 0.1 371 /journal.pone.0086756.g002 
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Figures 3. Box plots illustrating preoperative levels of urinary (A) norepinephrine and (B) epinephrine stratified accordingly to 
genotype clusters. 

doi:1 0.1 371 /journal.pone.0086756.g003 



Conclusion 

Somatic mutations are frequent events in PCC and PGL 
tumours. In patients with somatic carrier status there were no 
cases of recurrent nor metastatic disease. These findings suggest 
that analysis of tumour DNA could have an impact on the 
management of PCC and PGL patients and should be further 
investigated as prognostic factors in these diseases. 
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